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Nb can improve the resistance of Ni-based Hastelloy N alloy to Te-induced intergranular embrittlement. 
First-principles calculations are performed to research this mechanism by simulating the Ni(111) surface and 
the X` 5(012) grain boundary. The calculated adsorption energy suggests that Te atoms prefer diffusing along 
the grain boundary to forming the surface-reaction layer with Nb on surface of the Ni alloy. First-principles ten- 
sile tests show that the Nb segregation can enhance the cohesion of grain boundary. The strong Nb-Ni bonding 
can prevent the Te migration into the inside of the alloy. According to the Rice-Wang model, the strengthen- 
ing/embrittling energies of Nb and Te are calculated, along with their mechanical and chemical components. 
The chemical bonds and electronic structures are analyzed to uncover the physical origin of the different effects 
of Te and Nb. Our work sheds lights on the effect of Nb additive on the Te-induced intergranular embrittlement 


in Hastelloy N alloy on the atomic and electronic level. 
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I. INTRODUCTION 


Molten Salt Reactor (MSR) is the only liquid-fueled re- 
actor in the six most promising Generation IV reactor con- 
cepts [1]. As the structural material developed specially for 
MSR, Hastelloy N, a Ni-based alloy, has excellent corrosion 
resistance against molten salt, and was used in Molten Salt 
Reactor Experiment (MSRE) of the Oak Ridge National Lab- 
oratory (ORNL, USA). However, MSRE revealed that the 
usefulness of Hastelloy N is limited by its susceptibility to 
stress corrosion cracking (SCC) induced by Te, which is a 
most dangerous problem of Hastelloy N [2, 3]. Te, a fission 
product in fuel salt, tends to diffuse along the surface grain 
boundaries (GBs) of Hastelloy N and causes intergranular 
cracking eventually, which is related closely with SCC. 

To tackle this problem, a straightforward approach is to 
modify Hastelloy N by added alloying materials. MSRE 
found that adding Nb (1%~2%) to Hastelloy N was beneficial 
in reducing intergranular Te cracking, but still, it embrit- 
tled [4, 5]. The mechanism of the Nb effect on the Te- 
induced SCC is unknown: Nb may form a stable and innocu- 
ous telluride compound, or Nb hypothetically forms surface- 
reaction layers with Te in preference to the Te diffusion into 
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the alloy along the GBs, and so on [4, 5]. So, studies on this 
mechanism shall be helpful for developing more advanced Ni 
alloys, with adequate resistance to Te, for MSR. 

First-principles calculation is suitable to mechanism inves- 
tigations at atomic level. It was used successfully in study- 
ing effects of dopants or impurities in GB [6-9]. In this pa- 
per, we perform a first-principles calculation to clarify this 
mechanism by simulating a X` 5(012) Ni GB [7, 10] and the 
Ni(111) surface with the coexistence of Te and Nb. The re- 
sults about the effects of Te on Ni GB are in accordance with 
our previous work [9]. 


II. COMPUTATIONAL DETAILS 


Figure 1 is a schematic diagram of a X 5(012) Ni GB unit 
cell. It contains two reversely oriented grains with 80 Ni 
atoms. The atom layer is distinguished by the distance be- 
tween the layer and the GB plane. The GBO layer represents 
the hollow sites. There are four equivalent atomic sites in 
each layer. The geometry optimization calculation of GB was 
performed including cell optimization. We adopted the GB 
model in Ref. [7], and performed further optimization (in- 
cluding cell optimization) to find a more accurate GB model. 

The Ni(111) surface is modeled by a slab with a (4 x 4) 
surface periodic cell, which contains six layers of Ni atoms. 
The calculated lattice constant of bulk Ni used to build the 
Ni(111) surface is 3.52 A, which is in good agreement with 
the experimental result. The bottom layer without Te or Nb is 
fixed to its optimal bulk position to mimic the bulk. The vac- 
uum layer is about 12 A thick. The Nb-Ni(111) surface, with 
a Nb atom substituted for a Ni atom in the topmost layer of 
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Fig. 1. (Color online) Schematic diagram of (a) unit cell of 
>> 5(012) Ni GB. The atomic sites are labeled by numbers counted 
from the GB plane. For clarity, the gray and black balls represent 
atoms in layers with x = 0 (in the paper plane) and x = 0.25 (be- 
neath the paper plane) along the (100) direction, respectively. The 
other atoms with x = 0.5 and x = 0.75 are not shown. The right 
parts show the adsorption sites on the (111) surfaces of pure Ni (b) 
and Ni-Nb system (c): (1) (5) (1) top site; (2) bridge site; (3) hcp 
site; (4) fcc site. 


the Ni(111) surface (Fig. 1(c)), are also calculated to compare 
with the Ni(111) surface. 

Spin-polarized electronic state calculations were per- 
formed within the DFT [11, 12] using Vienna ab-initio sim- 
ulation package (VASP) [13]. Projector-augmented plane- 
wave (PAW) [14] methods were employed with the PBE gen- 
eralized gradient approximation (GGA) [15]. The wave func- 
tions were expanded in a plane-wave basis set with a cutoff 
energy of 350eV. The Brillouin zone was sampled using a 
3 x 3 x 1 k-point mesh. 

First-principles tensile tests were carried out to study the 
GB strength with Te or Nb in GB region. To simplify the cal- 
culations, the lattice dimensions in the GB plane were fixed 
to neglect the Poisson’s ratio. A uniaxial tensile strain was 
exerted in the GB normal direction (i.e. the (012) direction). 
In each strain step, the starting atomic configuration is taken 
from the relaxed configuration of the preceding step by an 
increment of 2% to ensure the continuous strain path. 


Il. RESULTS AND DISCUSSION 
A. Adsorption energy of Te 


The adsorption energies of Te, Faq, on the Ni(111) and Nb- 
Ni(111) surfaces are calculated by: 


Eaa = ETe-sub = E atom, Te = Egib (1) 
where Fte-sub, Fatom,Te and Esub refer to the calculated total 


energies of the optimized substrate with the adsorbate (i.e. 
a Te atom), one isolated Te atom, and the clean substrate, 
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respectively. A strongly negative value of Faq means intense 
binding between the Te atom and the substrate. 

As shown in Fig. 1, the top, bridge, hexagonal close-packed 
(hcp), and face-centered cubic (fcc) sites were considered for 
the adsorption. Table 1 shows the calculated adsorption en- 
ergies (in eV) of Te at each site on surface of Ni(111) and 
Nb-Ni(111), and the two data groups in three site types are 
roughly the same. These kinds of substitution of Nb for Ni do 
not change too much of the adsorption energy of Te, whereas 
the hcp site of Nb-Ni(111) is unstable for Te adsorption due 
to the existence of the Nb atom, and the Te atom is repelled 
to move from the hcp site to a farther site (Site 5 in Fig. 1(c)). 
As a result, Te atoms do not preferentially form the strong 
binding with Nb atoms on the surface of Ni-Nb alloy. So, 
the resistance of Nb to the Te-induced SCC in Hastelloy N 
cannot be attributed to the hypothetical formation of surface- 
reaction layers between Te and Nb. And Te would prefer to 
diffuse into the alloy along the GBs. The Nb effect in GB 
with the coexistence of Te will be discussed later. 


TABLE 1. Adsorption energies (eV) of Te at different sites on the 
Ni(111) and Nb-Ni(111) surfaces (See Fig. 1 for the atomic sites) 


Surface top bridge hep fcc 
Ni(111) —3.37 —4.13 —4.18 —4.26 
Ni(111)-Nb —3.12 —4.05 unstable —4.43 


B. First-principles tensile tests 


To further understand how the Nb additive affects the Ni 
GB at presence of Te, first-principles tensile tests were carried 
out to investigate the maximum strength of the GB and its 
fracture process. Te and Nb atoms, which are greater in size 
than Ni, prefer to occupy substitution sites (Site 1 in Fig. 1) 
rather than interstitial sites (Site 0 in Fig. 1) on the GB plane. 
So only the substitution case is considered. There are four 
sites in Layer 1. For simplification, the comparison was made 
among the clean GB, the GB + Nb layer (4 Nb atoms in Layer 
1), the GB + mixed layer of Nb and Te (2 Nb and 2 Te in 
Layer 1), and the GB + Te layer. 

As shown in Fig. 2, the GB + Nb layer has the largest ten- 
sile strength (21.3 GPa) at the strain of 28%, while that of 
the clean GB case is a little lower. On the other hand, the 
maximum strength of the GB + Te layer is about one-half of 
the case of GB + Nb layer. However, when the four sites 
in Layer | are occupied by 2 Nb and 2 Te atoms, the maxi- 
mum GB strength increases to 16.4 GPa at the strain of 20%, 
which is obviously improved compared with the GB + Te 
layer. Clearly, Nb segregation enhances the Ni GB cohesion, 
and Te in the GB region induces the Ni GB embrittlement, 
hence the inhibition of the Te-induced SCC in Ni GB by the 
segregated Nb atoms. 

In the region of strain < 28%, the elastic deformation oc- 
curs for the GB + Nb layer. From strain = 32%~36%, this 
GB undergoes the plastic deformation, and positions of the 
atoms in the GB cell are not layer by layer any longer. For 
the clean GB, GB + Te layer and GB + 2(Nb + Te), the frac- 
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Fig. 2. (Color online) Calculated tensile stress as a function of strain 
for the clean GB, GB + Nb layer (4 Nb atoms in Layer 1), GB + 
mixed layer of Nb and Te (2 Nb and 2 Te atoms in Layer 1), and GB 
+ Te layer. 


ture surface of the GB is indeed the GB plane. However, the 
fracture surface of the GB + Nb layer is not the GB plane but 
the plane between the Layers 2 and 3 (See Fig. 1). As a re- 
sult, the Nb(1)—Ni(2) bonds are stronger than the Ni(2)-Ni(3) 
bonds in the GB + Nb layer, and the corresponding Ni(1)- 
Ni(2) bonds in the clean GB. 


C. Strengthening/embrittling energy 


According to the Rice-Wang model [16], effects of the var- 
ious elements on the GB cohesion can be determined by the 
strengthening/embrittling energy, AF, which is defined as 


AE = (Ecp,doped — Eag) — (Ers.doped — Ers), (D) 


where, Egp.doped> EGB, EFs,dopea and Eps represent the total 
energies of the doped-GB, clean GB, doped-free surface (FS) 
and clean FS, respectively. A positive value of AE means 
embrittlement of the GB, and a negative value indicates en- 
hancement of the GB. 

To gain a deeper understanding, the strengthen- 
ing/embrittling energy can be decomposed into the me- 
chanical and chemical components. The procedure of 
decomposition in Refs. [10, 17—19] was used to perform the 
analysis. 

The calculated AF and its mechanical and chemical com- 
ponents for Te and Nb are listed in Table 2. According to the 
calculated strengthening/embrittling energy, Te, with a posi- 
tive value, is an embrittler, and Nb, with a negative value, is a 
cohesion enhancer. The mechanical components of Te and Nb 
are both positive. This is due to that the bigger atomic sizes 
of Te and Nb cause the GB expansion. On the other hand, the 
chemical component of Nb is strongly negative and plays a 
dominant role in the strengthening/embrittling energy. How- 
ever, Te has a small value of the chemical component, which 
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TABLE 2. Adsorption energies (eV) of Te at different sites on the 
Ni(111) and Nb-Ni(111) surfaces. (See Fig. 1 for the atomic sites.) 


Atoms AE Chemical Mechanical 
(eV) component (eV) component (eV) 

Te (this work) 2.16 0.14 2.01 

Te (Ref. [10]) 1.6 —0.2 1.8 

Nb (this work) —1.03 —2.89 1.86 


Nb (Ref. [20]) 


—1.05 - - 
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(b) GB with an Nb (Te) layer 


Fig. 3. (Color online) Calculated interatomic distances (A ) in the 
GB region for (a) clean GB, and (b) GB + 4 Nb in Layer 1 (a Nb 
layer) and the GB + 4 Te in Layer 1 (a Te layer). The data for the 
case of Te are given in parenthesis. The atoms are marked by their 
site numbers. 


has little effect on the strengthening/embrittling energy. So, 
the contrary effects of Te and Nb are mainly attributed to the 
difference between their chemical components. The results 
of our work agree well with the previous calculations [10, 20] 
(Table 2). 


D. Chemical bonds and electronic structures 


The atomic and electronic structures were studied for 
mechanisms of the different effects of Te and Nb on Ni GB. 
The calculated interatomic distances in the GB region are 
shown in Fig. 3 for the comparison. Being bigger than Ni 
in atomic radius, Te and Nb induce the GB expansion, as 
shown in Fig. 3. For example, comparing with in the clean 
GB, the Ni(3)-Ni(-3) distances in the GBs + Te layer and 
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(c) GB with an Nb layer 


Fig. 4. (Color online) Calculated charge density distribution (electron/Bohr?) in the (100) plane for (a) clean GB, (b) GB + Te layer and (c) 


GB + Nb layer. The atoms of interest are marked by their site numbers. 


GBs + Nb layer are elongated by 0.58 A and 0.24 A, respec- 
tively. The expansion can impair the GB cohesion, which 
is in accordance with the positive values of the strengthen- 
ing/embrittling energies for Te and Nb. Also, it can be seen in 
Fig. 3 that the Te-induced GB expansion is more serious than 
the Nb-induced expansion, as the larger positive mechanical 
component of Te is bigger than Nb. 

The chemical component of the strengthening/embrittling 
energy is thought to be induced by the charge redistribution 
due to the existence of the doped atoms [10]. With only a 
small difference in mechanical components for Te and Nb 
(see Table 2), their chemical components are comparable. 
Fig. 4 shows the calculated charge density distributions for 
the clean GB, GB + Te layer and GB + Nb layer. The Te(1)- 
Ni(2)/Ni(-2) bond in Fig. 4(b) and Nb(1)-Ni(2)/Ni(-2) bond in 
Fig. 4(c) are stronger than the corresponding Ni(1)-Ni(2)/Ni(- 
2) bond in Fig. 4(a), as judged by the charge densities along 
these bonds. But these strong Te-Ni bonds cannot enhance 
the GB cohesion that much, since directions of the bonds are 
almost in parallel with the GB plane. On the other hand, the 
Ni(2)-Ni(-2) and Ni(1)-Ni(4)/Ni(-4) bonds are normal to the 
GB plane and exert the main cohesive force to hold the two 
grains together as displayed in Fig. 4(a). Replacing Ni with 
Te in Fig. 4(b) makes the Ni(2)-Ni(-2) bond much weaker 
than the corresponding Ni(2)-Ni(-2) bond in Fig. 4(a). How- 
ever, for the GB with a Nb layer in Fig. 4(c), the Ni(2)-Ni(-2) 
bond is a little weaker than the corresponding Ni(2)-Ni(-2) 
bond in Fig. 4(a), but the Nb(1)-Ni(4)/Ni(-4) bond is obvi- 
ously stronger than the Ni(1)-Ni(4)/Ni(-4) bond in Fig. 4(a). 
As a result, the remarkable differences between the chemi- 
cal components of the strengthening/embrittling energies for 
Te and Nb are induced by these charge redistributions. The 
charge density for the GB + Te layer in Fig. 4(b) is in accor- 
dance with the result in Ref. [10]. 

As the conjugation interface between two misoriented 
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